BEC-BCS Crossover in the Nambu-Jona-Lasinio Model of QCD 
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The BEC-BCS crossover in QCD at finite baryon and isospin chemical potentials is investigated in 
the Nambu-Jona-Lasinio model. The diquark condensation in two color QCD and the pion conden- 
sation in real QCD would undergo a BEC-BCS crossover when the corresponding chemical potential 
increases. We determined the crossover chemical potential as well as the BEC and BCS regions. 
The crossover is not triggered by increasing the strength of attractive interaction among quarks 
but driven by changing the charge density. The chiral symmetry restoration at finite temperature 
and density plays an important role in the BEC-BCS crossover. For real QCD, strong couplings in 
diquark and vector meson channels can induce a diquark BEC-BCS crossover in color superconduc- 
tor, and in the BEC region the chromomagnetic instability is fully cured and the ground state is a 
uniform phase. 
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I. INTRODUCTION 

There exists a rich phase structure of Quantum Chro- 
modynamics (QCD), for instance, the deconfinement pro- 
cess from hadron gas to quark-gluon plasma and the tran- 
sition from chiral symmetry breaking to the symmetry 
restoration [l| at high temperature and/or baryon den- 
sity, the color superconductivity at low temperature but 
high baryon density 0, Q, and the pion superfluidity at 
low temperature but high isospin density [4j. The physi- 
cal motivation to study the QCD phase diagram is closely 
related to the investigation of the early universe, compact 
stars and relativistic heavy ion collisions. 

While the perturbation theory of QCD can well de- 
scribe the properties of the new phases at high temper- 
ature and/or high density, the study on the phase tran- 
sitions themselves at moderate temperature and density 
depends on lattice QCD calculation and effective models 
with QCD symmetries. While there is not yet precise 
lattice result for real QCD at finite baryon density due 
to the fermion sign problem Q, it is in principle no prob- 
lem to do lattice simulation for two color QCD at finite 
baryon density [g, Ilk la |9( and real QCD at finite isospin 
density [Tol. [TTL fl2| . The QCD phase transitions at finite 
baryon and/or isospin chemical potential are also inves- 
tigated in many low energy effective models, such as chi- 
ral perturbation theory [l^Ejjli OJ, 03 OB EI HI, 
linear sigma 
modellZ 
model[l3 



model 21 



Ml 



23fl, Nambu-Jona-Lasinio 



|2& HI S31iOprS IHI, random matrix 
33j . ladder QCD[34|. strong coupling lattice 
QCD[35l| and global color model [3r3|. 

It is generally expected that there would ex- 
ist a crossover from Bose-Einstein condensation(BEC) 
to Bardeen- Cooper- Shriffer condensation(BCS) for di- 
quarks at finite baryon density and pions at finite isospin 
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density. In the chemical potential region above but close 
to the critical value fi c for diquark condensate in two 
color QCD or pion condensate in real QCD, since the 
deconfinement does not yet happen, the system should 
be in the BEC state of diquarks or pions. On the other 
hand, at a sufficiently high chemical potential fi » /i c , 
the ground state of the system becomes a BCS superfluid 
where quark-quark or quark-antiquark cooper pairs are 
condensed. Therefore, there should be a crossover from 
BEC state to BCS state when the chemical potential or 
density increases. The BEC-BCS crossover, which is a 
hot topic in and beyond condensed matter and ultracold 
fermion gas[U Ell pi El EE E3, E3, E3 , was recently 
extended to relativistic fermion superfluid [H EE E3- 
In these studies, the BEC-BCS crossover is induced 
by increasing the attractive coupling among fermions. 
However, this is not the only way to trigger BEC-BCS 
crossover. In this paper, we will study the BEC-BCS 
crossover induced by density effect and chiral symmetry 
restoration. This density induced BEC-BCS crossover is 
like the one found in nuclear matter [48| and non-dilute 
fermi gas(49[. 

The features of BEC-BCS crossover have been widely 
discussed in condensed matter and ultracold fermion 
gas[H EB EB EO, EE El • We list here only the essential 
characteristics: 

1) There exist two important temperatures, one is the 
phase transition temperature T c for the superfluid at 
which the order parameter vanishes, and the other is the 
molecule dissociation temperature T*. In the BCS limit, 
there exists no stable molecule, and T* approaches to 
T c . In the BEC limit, all fermions form stable difermion 
molecules, and T* becomes much larger than T c . 

2) The chemical potential is equal to the Fermi energy of 
non-interacting fermion gas in the BCS limit but becomes 
negative in the BEC region. In the BEC limit, the ab- 
solute value of the chemical potential tends to be half of 
the molecular binding energy. This is true both near the 
critical temperature and in the superfluid ground state. 

3) The fermion momentum distribution changes signifi- 
cantly when we go from the BCS to BEC. In the BCS 
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limit, the distribution near the Fermi surface is still very 
sharp and similar to the one in non-interacting fermion 
gas. However, it becomes very smooth in the whole mo- 
mentum space in the BEC limit. 

To be specific, we express the typical dispersion of 
fermion excitations in a superfluid as 



E P = \ /(£p-^) 2 +A 2 , 



(1) 



where A is the superfluid order parameter. In non- 
relativistic homogeneous system, we have £ p = p 2 /(2m). 
In the BCS case, we have /x > 0, the minimum of the 
dispersion is located at nonzero momentum |p| = y/2mfi 
and the excitation gap is A. However, when \i becomes 
negative in the BEC region, the excitation gap becomes 
y/ fi 2 + A 2 rather than A itself, and the minimum of the 
dispersion is located at |p| = 0. This can be regarded as 
the definition of a BEC-BCS crossover at zero temper- 
ature. How can we extend this definition to relativistic 
systems? In relativistic case, we have £ p = \/p 2 + m 2 . 
In comparison with the non-relativistic result, we should 
define a new chemical potential fijy = fi—m. For /i > m, 
we have /xjv > 0, the minimum of the dispersion is located 
at nonzero momentum |p| = \//i 2 — m 2 and the excita- 
tion gap is A. On the other hand, for // < m, /ijv becomes 
negative, the excitation gap becomes y/ JA^ + A 2 and the 
minimum of the dispersion is located at |p| = 0. This 
analysis indicates that the BEC-BCS crossover in rela- 
tivistic fermion gas is controlled by /xjv = (J-—m rather the 
chemical potential /i itself. An interesting phenomenon 
then arises: The BEC-BCS crossover would happen when 
the fermion mass m varies in the process of chiral sym- 
metry restoration at finite temperature and density. 

The effective models at the hadron level can only de- 
scribe the BEC state of hadrons and diquarks [Hoi [5l| . 
they can not describe the possible BEC-BCS crossover 
when the chemical potential increases. One of the ef- 
fective models that enables us to describe both quark 
and meson properties is the Nambu-Jona-Lasinio(NJL) 

models applied to quarks[El iHl iHE IHE IH3- Even 
though there is no confinement in the NJL model, the 
chiral phase transition line [53, |54| . l55l . I56I . [57], 15a . [59( in 
the temperature and baryon chemical potential (T — hb) 
plane calculated in the model is very close to the one ob- 
tained from lattice QCD. It is natural to extend the NJL 
model to studing diquark condensation at finite baryon 
chemical potential and pion condensation at finite isospin 
chemical potential. While there is no reliable lattice re- 
sult for real QCD with diquark condensation at finite 
baryon chemical potential, the calculation of diquark con- 
densation in the NJL model [28| agrees well with the lat- 
tice simulation of two color QCD0, 0. Also, the NJL 
model calculation of pion superfluidity [3(| agrees well 
with the lattice simulation of real QCD at finite isospin 
chemical potential HH, EH- ^ is natural to ask: How 
can such a model with quarks as elementary blocks de- 
scribe the BEC-BCS crossover of diquarks and pions? In 
non-relativistic fermion gas, the microscopic model with 



four-fermion interaction can describe well the BEC-BCS 
crossover at least at zero temperature [HI, Hi]]. Motivated 
by this fact, we believe that the NJL model should de- 
scribe well the BEC-BCS crossover in QCD at finite den- 
sity. 

The paper is organized as follows. In section [H] we in- 
vestigate the diquark BEC-BCS crossover at finite baryon 
chemical potential in the NJL model of two color QCD. 
The possible diquark BEC-BCS crossover in real QCD is 
investigated in the NJL model in section [TTT1 In section 
IIVI we study the BEC-BCS crossover of pion condensa- 
tion at finite isospin chemical potential. We summarize 
in section [V] 



II. DIQUARK BEC-BCS CROSSOVER IN TWO 
COLOR NJL MODEL 



Let us first consider the diquark condensation in two 
color QCD. The advantage to take two color QCD is that 
the diquarks in this case are colorless baryons and the di- 
quark condensation breaks the baryon symmetry Us(l) 
rather than the color symmetry SUc(2). The confine- 
ment in two color QCD is less important than in three 
color QCD. We start with the two color and two flavor 
NJL model 



(2) 



where mo is the current quark mass, and U are the 
Pauli matrices in flavor and color spaces, and the two 
coupling constants G s and Gd are connected by a Fierz 
transformation in color spacefill], G s — Gd — G. In 
the chiral limit, the Lagrangian ^ with flavor num- 
ber Nf = 2 has an enlarged flavor symmetry SU(2iVf), 
which is called Pauli- Guersey symmetry [28| and connects 
quarks and anti-quarks. As a consequence, diquarks are 
color singlet baryons and diquarks and pions become de- 
generate. 

The key quantity describing a thermodynamic system 
is the partition function Z which can be expressed as 



Z 



[dip][dip]e 5 ° ' dr J '<* 3 x(£+£f-ih<oV) ( 3 ) 



in the imaginary time formulism of finite temperature 
field theory, where the baryon chemical potential /is is in- 
troduced explicitly, and [3 is the inverse temperature, (3 = 
l/T. Using the Hubbard-Stratonovich transformation, 
we introduce the auxiliary meson fields a = —2G^>^>, it = 
~2G^ji r y5Ttp and diquark field <f> = —2Gip T Cij5T 2 t 2 tp, 
and the partition function can be written as 



Z = [dV] [d&] [da] [dn] [defy] [d<j>* ] 



X rfr/d 3 x£ cff 



with the effective Lagrangian 



(4) 



(5) 
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where we have introduced the Nambu-Gorkov spinors 



and the kernel K, is defined as 



, ¥ = ( 4, ) , (6) 



/C[(T, 7T, ( 



M+ 175^2*2 



(7) 



with M± — i^dfj, — mo ± /IJ370/2 — (cr ± i7sr • 7r). Inte- 
grating out the quark degrees of freedom, we obtain 



Z = / [da}[dTt}[d(f>}[d(l)*}e 



(8) 



with the bosonic effective action 
Jo 



Stff [c, 7T, . 



4G 



1 



Trln/C[cr, n, < 



(9) 



A. Diquark Fluctuation at T > T c 

The diquarks would be condensed when the baryon 
chemical potential /173 is larger than their mass = 
m w . In this subsection we focus on the region above 
the critical temperature T c where the diquark condensate 
vanishes. 

After a field shift a —> (a)+a with (cr) = -2G(^),the 
effective action S^g at zeroth order in meson and diquark 
fields gives the mean field thermodynamic potential f2, 

n 1 c(o) (m-m ) 2 1 ! 
n= ^ = 4G + W lndet5 ' (10) 

where m — m + (a) is the effective quark mass, and S 
is the quark propagator at mean field level 



S = (V<9p - m + ^r7o°- 3 ) 



(11) 



with <Ti being the Pauli matrices in the Nambu-Gorkov 
space. The quadratic term of the effective action reads 



4G 



-Tr{«SS[cr,7r,</»]«SS[cr,7r, (12) 



where £ is defined as 



£[<7, 7T, 1 



(13) 



«750*T 2 i2 cr - fy s r • 7r 
In momentum space, the second order effective action 



can be expressed as 



2G 



2G-^ (fc) 

^-n-(*) 



2G 



n ff (fc) cr(-A;)cr(A;) 

7T(-fc) • 7T(fc) 

<t?{-k)(t>{k) 

<K-k)4T{k) , (14) 



where fc = (fc ,k) with fc = iu> n — 2imtT(n = 
0, ±1, ±2, • • • ) is the meson or diquark four momentum, 
and J2k ~ "iTj^nJ rf 3 k/(27r) 3 indicates integration over 
the three momentum k and summation over the fre- 
quency Lo n . The polarization functions for the mesons, 
diquark and anti-diquark can be evaluated as 

U a {k) = 2IIi(fc; mb), IMA;) = 2H 2 (k; Ms ), 
U d (k) = 2n 3 (fc; y. B ), n d -(fc) = 2n 4 (fc; Ms) (15) 

with the functions 111^^4 (fc; /_i) listed in Appendix A. 

Without loss of generality, we consider the case [Ib > 
0. The transition temperature T c for the diquark conden- 
sation is determined by the well-known Thouless criterion 



1 - 2Gn d (fc = 0) 







T=T C 



(16) 



together with the gap equation for the effective quark 
mass derived from the first order derivative of the ther- 
modynamic potential, 



mo 



16Gm 



d 3 P l-f(E+)-f(E-) 



(2tt) 3 



(17) 



T=T C 



with the particle energies E p = \/p 2 + m? and E p = 
E p ±fi B /2. 

The NJL model is non-renormalizable, a proper regu- 
larization is needed to avoid ultraviolet divergence. For 
instance, one can add a form factor in the above momen- 
tum integrations. Since our goal is to study the BEC- 
BCS crossover which happens far away from the asymp- 
totic region, we employ, for the sake of simplicity, a hard 
three momentum cutoff A to regularize the above equa- 
tions. In the following numerical calculations, we take 
the current quark mass mo = 5 MeV, the coupling con- 
stant G = 1.5 x 4.93 GeV~ 2 and the cutoff A = 653 
MeV to fit the pion mass = 134 MeV, the pion decay 
constant f % — 93 MeV and the constituent quark mass 
to = 300 MeV in the vacuum. 

The critical temperature as a function of baryon chem- 
ical potential is shown in Fig[TJ The diquark condensed 
phase starts at fis — and the critical temperature 
increases with [Ib until a saturation value which is about 
the critical temperature for chiral symmetry restoration 
at [Lb = 0, Tq = 185 MeV. At moderate baryon chemical 
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potential, the critical temperature can be well described 

El 



T C = T \ 1 



77V 
MB 



(18) 
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FIG. 1: The superfluid phase transition temperature T c scaled 
by To as a function of baryon chemical potential /is scaled by 

TO.. 



We discuss now the mesons and diquarks above T c . 
The energy dispersions of the mesons and diquarks are 
defined by the poles of their propagators, 



l-2GIIi(w(k),k) =0, 



, d, d. 



(19) 



While the full dispersion laws can be obtained numeri- 
cally, we are interested here only in the meson and di- 
quark masses. For a and it, they do not carry baryon 
number, and their masses are defined as m a — ^(O) 
and = w ff (0). However, diquarks and anti-diquarks 
carry baryon numbers and their masses are not exactly 
their dispersions at k = 0. To define a proper di- 
quark or anti-diquark mass, we should subtract the corre- 
sponding baryon chemical potential from the dispersion, 
m d = Lu d (0) + [i B and m d = wj(0) - Ms- For \i B > 0, by 
comparing equation (|16p for the transition temperature 
T c with the above definition of diquark mass, the diquark 
mass at T c is equal to the baryon chemical potential, 



m d (T c ) = fx B , 



(20) 



which is consistent with the physical picture of relativistic 
BEC in boson field theory. In Fig[2]we show the meson 
and diquark mass spectrum above the critical tempera- 
ture for fiB < rn-x and fin > Tn w . For fiB < Tn w , the 
meson mass spectrum is similar to the case at \j,b = 0, 
but the diquark mass is slightly different from the anti- 
diquark mass at high temperature. For \xb > m^, the 
mass spectrum starts at T c where the diquark mass is ex- 
actly equal to [Ib- For any we found that the meson 
or diquark mass becomes divergent at a limit temper- 
ature T* which indicates the dissociation of the meson 
or diquark resonances. For small baryon chemical po- 
tential, this dissociation temperature is about two times 



the critical temperature for chiral symmetry restoration, 
which is consistent with the results in 61]. With in- 
creasing baryon chemical potential, while the dissocia- 
tion temperature for diquarks and mesons decreases, the 
one for anti-diquarks increases. Since the diquarks are 
condensed at positive /j,b, the baryon chemical potential 
dependence of the diquark dissociation temperature indi- 
cates a crossover from diquark BEC to BCS superfluidity. 
In Fig[3l we show the superfluid transition temperature 
T c and the dissociation temperature T* for mesons, di- 
quarks and anti-diquarks. Any T* is much larger than T c 
at small /ib, but for mesons and diquarks the two tem- 
peratures tend to coincide with each other at sufficiently 
large \x B . 
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FIG. 2: The meson and diquark masses as functions of tem- 
perature for fiB = 0.06 GeV and 0.2 GeV. 



0.5 



0.4 



0.3 



0.2 



0.1 



T* 






-• r .l~- ■ 


. 1 . i G,7r 








" " * * 








c 






) 0.2 


0.4 0.6 


0.8 




H B (GeV) 





FIG. 3: The meson and diquark dissociation temperatures as 
functions of the baryon chemical potential. 
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To investigate the diquark excitation in detail, we con- 
sider the diquark spectral function 



p(w,k) = -2IiriD fl (w,k), 



(21) 



where Dfl(w,k) = T>(uj + imk) is the analytical contin- 
uation of the diquark Green function 2?(iu; n ,k) defined 
as 



r>(iw„,k) 



2G 



(22) 



1 - 2GU d (tuj n ,k) 
The spectral function p can be rewritten as 

p(w,k)= (23) 

-8G 2 Imn d (w -Hrj,k) 
[1 - 2GReU d {uj + irj, k)] 2 + [2GTmII d (u) + ir), k)] 2 ' 

At zero momentum k = 0, the real and imaginary parts 
of the diquark polarization function can be evaluated as 



Ren d (w -I- i77,0) = 



(24) 



d 3 p 



l-2/(£-) l-2/(^ 



(2tt) 3 
ImfldttJ + ir\, 0) = 
d 3 p 



uj — 2E P + (is w + 2E P + pb 



-8tt 



(2tt) 3 



{l-2.f{E v ))5{Lo-2E v + l i B ) 
-(l-2f(E+))5(cj + 2E p + fi B ) 



- [ Pu} E Pul (1 - 2f(E~J)] 6(w + /is — 2m), 



with 



Pa 



— m . 



(25) 



For iib < 2m, a diquark can decay into two quarks only 
when its energy satisfies uj > 2m — ilb, due to the step 
function O in the imaginary part. However, for > 2m, 
a diquark can decay into two quarks at any energy and 
hence becomes an unstable resonance. In FigUwe show 
the diquark spectral function at zero momentum for sev- 
eral values of baryon chemical potential near the criti- 
cal temperature T c . For small baryon chemical potential 
fiB < 2m, the diquarks are in stable bound state at small 
energy u> < 2m — \xb and in unstable resonant state at 
large uj > 2m — iib- For large baryon chemical poten- 
tial hb > 2m, the diquarks are impossible to stay in 
bound state. With increasing baryon chemical potential, 
the diquarks become more and more unstable and finally 
disappear. This naturally supports the picture of diquark 
BEC-BCS crossover. 

To understand more clearly the crossover from di- 
quark BEC to BCS superfluidity, we come back to the 
BEC-BCS crossover in non-relativistic condensed matter 
physics. A signal of the crossover is that in the BEC 
region the chemical potential becomes negative and its 
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FIG. 4: The scaled diquark spectral function at zero momen- 
tum for small and large baryon chemical potentials. 
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FIG. 5: The effective quark mass and the corresponding chem- 
ical potential ^ as functions of baryon chemical potential at 
the critical temperature. 



absolute value tends to be half of the molecule binding 
energy. Therefore, the fermions are heavy and hard to 
be excited even at finite temperature, and only molecules 
survive and condense in the BEC region. In relativistic 
systems, we should subtract the fermion mass from the 
chemical potential in order to compare with the non- 
relativistic result. For this purpose, we define a new 
chemical potential 



Mb 

ri N = — m. 



(26) 



In the deep BCS limit with /i^ — > oo, we have m — > and 
^ Mb/2, but in the deep BEC limit with \xb — > m„, 
m tends to be its vacuum value, and /xjv becomes neg- 



6 



ative and its absolute value approaches to half of the 
diquark binding energy — 2to — m,. In FigfS] we show 
the effective quark mass m and the chemical potential 
/ijv at the critical temperature. The effective quark mass 
decreases gradually from its vacuum value to zero with 
increasing baryon chemical potential. Since the lowest 
quark excitation energy is -\/p 2 + to 2 — /is/2, the quarks 
are hard to be excited near the critical temperature if to 
is much larger than ^b/2. Obviously, at small /ib, the 
diquarks are easy to be excited near T c . This means that 
even though the elementary blocks of the NJL model are 
quarks, the true physical picture at small baryon chemi- 
cal potential is the BEC of diquarks. The baryon chem- 
ical potential fj, B corresponding to the crossover can be 
defined as [i B — 2to which is just the point where the 
diquarks become unstable resonances. From FigO /i B is 
about 240 MeV. We will give an analytical expression for 
/i B in the next subsection. 



B. Diquark Condensation at T < T c 

In this section we focus on the quark and diquark be- 
havior in the superfluid phase. We start from the effec- 
tive action Since diquarks are condensed, we should 
introduce not only the chiral condensate (a) but also the 
diquark condensates, 



(27) 



A nonzero diquark condensate A ^ means spontaneous 
breaking of Us(l) baryon number symmetry in the two 
color QCD. The phase factor 8 indicates the direction 
of the Us(l) symmetry breaking. For a homogeneous 
superfluid, we can choose 9 — without loss of general- 
ity. A gapless Goldstone boson will appear, which can 
be identified as the quantum fluctuation in the phase di- 
rection. This Goldstone boson is just the Bogoliubov 
phonon in the superfluid. 

After a field shift a — > (a) + a and <f> — > A + <j>, the 
effective action at zeroth order in meson and diquark 
fields gives the mean field thermodynamic potential 



o = _L c(°) = (m-m ) 2 + A 2 
PV cff 4G 



27^1ndetS A \ 

(28) 

where 5a is the mean field quark propagator in the su- 
perfluid phase. In momentum space it can be evaluated 
as 

Sa(p) = (VJV ~m+ ^y7o^3 + H 5 AT 2 t 2 a^j ' (29) 

with the quark four momentum p = (iv n ,p) = ((2n + 
l)wrT,p) (n = 0,±f,±2, • • •)• In the Nambu-Gorkov 
space, it can be explicitly expressed as a matrix 



Sa(p) 



Sn{p) S 12 (p) 
Snip) S 2 2(p) 



(30) 



with the elements 

(iv n + E~) A+7o (iv n - E+) A_ 7o 



S 22 
S\ 2 
S 2 \ 



{iu n ) 2 - {E- A ) 2 



{iv n - E p ) A_7q (iv n + E+) A +7o 

(«/„)* -(££) 2 + 
-iAT 2 i 2 A + 7 5 



+ \2 



-^At 2 £ 2 A_ 75 

{iv n y-{E^)^ {iv n y-(Etr 

-iAr 2 t 2 A_75 -iAr 2 t 2 A +75 



(«"») - 

where A± are the energy projectors 



(iv n y-(E+y- 



A ± (P) = \ 



I ± 



7o (7 - P + m ) 
E n 



(31) 



(32) 



and E 



(£p) 2 



A 2 are quark energies in the su- 



perfluid phase. The momentum distributions of quarks 
and anti-quarks can be calculated from the positive and 
negative energy components of the diagonal propagators 
S11 and S 22 , 



n q (p) 



d,, ,2 



iv n + E' 



««(p) = 



- {iv n f - (E A ) 
x - iv n + E. 



1 f E~ 

n*nW — [ 1 — 

2 I E 



A 



1 ( E i" 

^(^-(^ el ^ = H 1 ~^ ) - ,:>,:!! 



It has been demonstrated in non-relativistic models 
that the BCS mean field theory can describe well the 
BEC-BCS crossover at low temperature, T <C T c . Here 
we will treat the ground state in the mean field approx- 
imation. The gap equations to determine the effective 
quark mass to and diquark condensate A can be obtained 
by the minimum of the thermodynamic potential, 



= 0, 



= 0. 



(34) 



dm ' OA 
With the explicit form of the thermodynamic potential 



= 



(to - to ) 2 + A 2 
4G 



d 3 p 
(2^)3 



[((E+)+aE±)], 
(35) 



where £ is defined as ((x) = x/2 + (3 1 ln(l + e f3x ), we 
obtain the gap equations at zero temperature 



to — toq 



A 



8Gr 



8GA 



d 3 p 1 (E. 



(2n) 3 E p \E t 
d 3 p 



Ei 



(2tt) 3 \E 



El 



(36) 



Numerical solution of the gap equations is shown in FigJS] 
Our results for the chiral and diquark condensates agree 
quite well with the lattice data[28|. For [i B < the 
ground state is the same as the vacuum state and the 
baryon density keeps zero. For /is > m n the diquark 
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condensate and baryon density become nonzero. The 
critical baryon chemical potential /Xg is exactly the di- 
quark mass md in the vacuum, and the phase transition 
is of second order. Since the chiral symmetry is sponta- 
neously broken at small /ib and almost restored at large 
Hb , the effective quark mass m plays an important role at 
small fiB but can be neglected at large fiB- In fact, in the 
region above but close to the critical value \x c B , the chiral 
condensate, effective quark mass and the diquark con- 
densate as functions of /x B can be well described by [13] 



Ms 




<7>(0) 



Ms 



(37) 



which are consistent with the results from the chiral ef- 
fective theory fill]. 
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FIG. 6: The diquark condensate and the effective quark 
mass scaled by the quark mass in the vacuum as functions 
of \ib /m n at zero temperature. 



To explain the BEC-BCS crossover, we further con- 
sider the dispersions of quark excitations, 



VP 2 



±Mb/2 +A2 



(38) 



In the case of [Ib > 0, E^ is the anti-particle excita- 
tion. In Fig[7] we show the dispersion E^ at Mb = 
0.15, 0.5, 0.9 GeV. Obviously, the fermion excitations 
are always gapped. At small \xb with Mb/2 < m(fiB), the 
minimum of the dispersion is at |p| = where the energy 
gap is a/maj- + A 2 with the corresponding non-relativistic 
chemical potential fijy = Mb/2 - m introduced in the last 
subsection. However, at large [Ib with mb/2 > m(fiB), 
the minimum of the dispersion is shifted to |p| ~ Ms/2 
where the energy gap is A. Such a phenomenon is a 
signal of the BEC-BCS crossover. 

The BEC-BCS crossover is also reflected in the mo- 
mentum dependence of the quark occupation number. 
In Fig[7]we show also the quark momentum distribution 
n q (p) at [i B — 0.15, 0.5, 0.9 GeV. At small [i B such as 
[iB = 0.15 GeV the occupation number is very small and 



smooth in the whole momentum region, while at large fiB 
the occupation number becomes large near |p| = and 
drops down with increasing momentum rapidly. When 
Hb becomes very large, the occupation number is indeed 
of the BCS type. 




FIG. 7: The quark energy dispersion and momentum distri- 
bution at several baryon chemical potentials at zero temper- 
ature. 



The numerical results of the chemical potential /i^v is 
shown in FigO The turning point from negative jiN 
to positive fiN is at about fi B ~ 230 MeV, which is 
nearly the same as we obtained at the critical temper- 
ature shown in the last subsection. In fact, we can get 
an analytical expression for [i B using equation (|37p . The 
BEC region is obtained by requiring /in < 0, namely 



Ms 
2 



< m(/ii 



m(0) 



77V 

Ms 



from which we obtain 



Ms = 2m(0)m^ 



21 1/3 



(39) 



(40) 



it depends on the pion mass and effective quark mass in 
the vacuum. With the above chosen parameter set, we 
have ii b = 230 MeV. When the effective quark mass in 
the vacuum varies from 300 MeV to 500 MeV, fi B changes 
from 230 MeV to 270 MeV. Note that the explicit chi- 
ral symmetry breaking induced by nonzero current quark 
mass mo plays an important role in the study of BEC- 
BCS crossover. In the chiral limit with mo = 0, pions 
are massless and there will be no BEC state. 
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FIG. 8: The effective quark mass and non-relativistic chemical 
potential (in as functions of baryon chemical potential at zero 
temperature. 



We now turn to the meson and diquark excitations. 
The second order effective action which is quadratic in 
the meson and diquark fields can be expressed as 



dT / d J X- 



4G 



(41) 



~4 Tr { 5 A£[cr,7T,(/>]SA£[er,7r,0]} , 



and can be evaluated as 



2 ^ 

k 



2G 



-n«(fc) 



ipi{-k) Vj {k) (42) 



in momentum space, where ipi stand for the meson and 
diquark fields a, n, <f> and <p* , and the polarization func- 
tions ili, are defined as 



J ^Trp^Ab + ^r^^)] 



(43) 



with Ti being the interaction vertex defined in the La- 
grangian density @. 

Due to the off-diagonal elements of the quark propaga- 
tor (|5D|) . those eigen modes of the system above T c are, 
in general case, no longer the eigen modes of the Hamil- 
tonian in the superfluid phase. The new eigen modes 
are linear combinations of the old eigen modes and their 
masses are controlled by the determinant of the meson 
and diquark polarizations 



det 



2G 



iWfco = M,k = 0) 



0. 



(44) 



m and A coexist and are both large, but the mixing can 
be neglected at large baryon chemical potential where m 
approaches to zero. 



C. Phase Diagram in T — (ib Plane 

We can now summarize the above results and propose 
a phase diagram of two color QCD in the T — (Ib plane. 
The phase diagram is shown in FigJS] At low temperature 
and low baryon chemical potential, the matter should be 
in the normal hadron state. The thick dashed line at 
high temperature means the estimated phase transition 
from hadron gas to quark gas, which can not be calcu- 
lated in the NJL model but should exist in the system. 
When the baryon chemical potential becomes larger than 
the pion mass in the vacuum, the diquark BEC appears 
and exists up to another critical baryon chemical poten- 
tial (i B indicated by the vertical dashed line. At high 
enough baryon chemical potential, the matter will be- 
come the BCS superfluid where the quark-quark Cooper 
pairs are condensed. Between the BEC and BCS states 
there should exist a crossover region, like the pseudogap 
regime in high temperature superconductor and ultra- 
cold fcrmion gas. Note that the transition from diquark 
BEC to BCS superfluid is a smooth crossover. Since two 
color QCD can be successfully simulated on lattice, such 
a BEC-BCS crossover can be confirmed by measuring 
the quark energy gap and comparing it with the diquark 
condensate. The pseudogap phase at high temperature 
can also be confirmed by investigating the quark spectral 
function. 
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FIG. 9: The proposed phase diagram of diquark condensation 
in two color QCD in the T — (Ib plane. 



It can be analytically proven that while pions do not 
mix with the others and hence are still the eigen modes 
of the system, there is indeed a mixing among sigma, 
diquark and anti-diquark(28j, and this mixing leads to 
a gapless Goldstone boson. We can show that H a d and 
II^j are proportional to mA and n d j is proportional to 
A 2 . Therefore, the mixing between sigma and diquark 
or anti-diquark is very strong in the BEC region where 



III. DIQUARK BEC-BCS CROSSOVER IN 
THREE COLOR NJL MODEL 



Let us now consider the realistic world with color de- 
crees of freedom N c = 3. The three color NJL model 
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including the scalar diquark channel is defined as 



+G d 



a=2.5,7 



(^ 7 5 r 2 A a C^ T ) (^ T « 7 5 T 2 A a V) 



(45) 



where A a are the Gell-Mann matrices in color space. Dif- 
ferent from two color QCD, the scalar diquarks in three 
color case are color anti-triplets. Therefore, in real QCD, 
the problem of confinement becomes important. In this 
section, we will discuss the possible diquark BEC-BCS 
crossover without considering the effect of confinement. 



653 MeV. For convenience, we take the coupling ratio 
V = Gd/G s instead of G d . With the above parameter 
values, we find r) ma . x ~ 1.55 and rj m i n — 0.82. For other 
possible parameter values, the lower limit is roughly in 
the region 0.7 < rj min < 0.8. Physically speaking, if 
f] > ?7max, the vacuum becomes unstable, and if rj < rj m i n , 
a scalar diquark can decay into two quarks and becomes 
a unstable resonance. 

The t) dependence of the scalar diquark mass is shown 
in FigfTDl As estimated in Q, a scalar diquark is in a 
deeply bound state when the binding energy is in the re- 
gion 200 MeV < 2m — m d < 300 MeV, which corresponds 
to the diquark mass 300 MeV < m d < 400 MeV and to 
the coupling ratio 1.3 < rj < 1.4. 



A. Scalar Diquark Mass in the Vacuum 

We first discuss the scalar diquark mass in the vacuum 
with T = [Ib = 0. In the random phase approximation, 
the diquark mass is determined by the pole equation 

l-2G d n d (A;o=m d ,k = 0)=0, (46) 

where the diquark polarization function is denned as 

-^JLTr{i l5 S(p + k)i l5 S(p)} (47) 

with the mean field quark propagator S(p) — ("y^p^ — 
to) -1 . The effective quark mass m defined as m — mo — 
2G s ('ipip} satisfies the gap equation 

d 3 p 1 



m — too = 2AG S 



(2^) 3 ^ 



(48) 



Taking the trace in the Dirac space and summation 
over the quark frequency, the diquark mass mj is con- 
trolled by 



1 = 8G d 



1 



1 



(2tt) 3 V E p + m d /2 E p - m d /2 



(49) 



If diquarks can exist as stable bound states in the vac- 
uum, their mass must satisfy the constraint < m d < 
2m. As a consequence, the coupling constant in the scalar 
diquark channel should be in the region G™ m < G d < 
G™ ax [12, EH , where the upper limit is nearly model in- 
dependent, 



Gin; 
d 



Ig s -^ 

2 m — mp 



2 



(50) 



but the lower limit depends on the model parameters, 

^2 



Gmin 
,l — 



(Av/A 5 



to 2 + to 2 In A +^T^ 



(51) 



In the three color NJL model, the model parameters 
are set to be to = 5 MeV, G s = 4.93 GeV -2 and A = 
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FIG. 10: The scalar diquark mass rrid in the vacuum as a 
function of the coupling ratio rj — Gd/Gs- 



B. Diquark BEC-BCS crossover 

At large enough baryon chemical potential, the di- 
quarks will condense. The diquark condensate is defined 
as 



-2G d (ip T Cij 5 T 2 \ a ^}. 



(52) 



Due to the color SU(3) symmetry, we can choose a spe- 
cific gauge A 2 = A 7^ 0, A 5 = A 7 = 0. In this gauge, 
the red and green quarks participate in the condensa- 
tion, but the blue one does not. We firstly consider the 
simplest case where the chemical potentials for the red, 
green and blue quarks are equal, \x r = \x g = /if, = /k_b/3. 
General cases with color neutrality and vector meson cou- 
pling will be discussed in the following subsections. The 
thermodynamic potential at mean field level reads 64 



O = 



(to - to ) 5 



4G S 



A 2 

4G d 



(53) 



d 3 p 



2 (C(E+) + ((E A )) + C(E+) + C(£ 6 ")] , 



(2tt) 3 

where the quark energies are defined as 



(Ep ± /W3) + A 2 



E± = E p ±ji B /3. 



(54) 
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Minimizing the thermodynamic potential, we obtain the 
gap equations for the chiral and color condensates m and 
A at zero temperature, 



m — mo 



= 8G S 



A = 8G d A 



d 3 p 1 

(2tt)3 U; 



5 



(K) 



A 



1 

si 



(55) 



The above gap equations are almost the same as in 
the two color case, except the term Q(E b ') in the first 
equation. For i] < 7i m i n , diquark condensation takes place 
at about \x c B ~ 3m and the phase transition is of first 
order. In this case, there exists no BEC region where 
fiN = fi r — m = fi B /3 — m is negative. On the other hand, 
For ?y m i n < 77 < 77 max , the diquarks become condensed at 
fi L B = 3m d /2 and the phase transition is of second order. 
The proof is as follows. For [i B < 3m d /2 < 3m, the gap 
equation for m keeps the same form as in the vacuum, 
and the gap equation for A becomes 



l = 8G d 



1 



1 



d 3 i 



(2tt)3 \E p + f i B /3 £ P -M c B /3 



(56) 



at the phase transition of color superconductivity with 
A = 0. From the comparison with the diquark mass 
equation in the vacuum, the critical baryon chemical 
potential for color superconductivity should be fi c B = 
3m<j/2[63]. In this case, there must exist a BEC region 
where fj,N is negative, since at the phase transition point 
we have /in = ^ C B /'S — m = m d /2 — m < 0. 

In Fig [Til we show the effective quark mass and di- 
quark condensate as functions of baryon chemical po- 
tential. From the behavior of fi^, there exists indeed 
a diquark BEC region at intermediate baryon chemical 
potential. At sufficiently large baryon chemical potential 
the diquark condensate is in the BCS form. 

In the diquark BEC region, we assume that the effec- 
tive quark mass behaves as 



m(0) 



Mb 
/is 



a > 0. 



(57) 



From our numerical calculation, the value of a depends 
on the coupling 77. Taking m(/Js) = /jb/3 from /j^ = 
at the end point of the BEC, the diquark BEC region 
ends at 



Ms 



[2m(0)m^ 



Therefore, the interval of /iB for diquark BEC is 

'2m(0) 



A// B 



Ms 



Mb 



;m d 



m d 



1 



(58) 



(59) 



The largest BEC region takes place at md — 185 MeV 
corresponding to rj ~ 1.46. 

The critical temperature T c can be determined by solv- 
ing the gap equations with A = 0. Above the critical 
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FIG. 11: The diquark condensate A, effective quark mass m 
and the chemical potential as functions of baryon chemical 
potential for several values of the coupling ratio 77 at zero 
temperature. The color neutrality is not taken into account. 



temperature, all scalar diquarks become degenerate and 
their polarization function reads Tl d (k) = Il3(fc;2/^B/3) 
with II3 listed in Appendix |A"1 Similar to the calculation 
in the two color case, the diquarks are in bound state in 
the BEC region with /xb/3 < m, and become unstable 
resonances in the BCS region. The numerical results and 
discussions are quite similar to that in the two color NJL 
model. 



C. Effect of Color Neutrality 

In three color QCD, diquarks are no longer colorless 
and the requirement of color neutrality is not automat- 
ically satisfied once diquarks are included. This can be 
seen if we compute the expectation values of the color 
charges (Q a ) = (-tpXaJoip)- In the gauge with A 2 ^ 
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and A 5 = A7 = 0, we find (Qg) 7^ 0. To solve this 
problem, we can introduce a color chemical potential fig 
corresponding to the 8-th color charge. It may be dynam- 
ically generated by a gluon condensate (Ag). The color 
chemical potential enters the gap equations by replacing 
the quark chemical potentials fi r = fj, g = fib = ps/3 
by fi r = fig = ps/3 + Ms/3 and fib = ps/3 ~ 2p 8 /3. 
The color charge neutrality condition is guaranteed by 
ng = —dfl/dfis — 0, namely 



d 3 p 



Ei 



5 

ET. 



26 



(60) 



Taking into account the color neutrality, the recalcu- 
lated effective quark mass and diquark condensate are 
shown in Fig(T2] Comparing the behavior of the effective 
chemical potential /i/v = \i r — m in the two cases with 
and without considering color neutrality, we find that the 
requirement of color neutrality disfavors diquark BEC. 
However, it does not cancel the BEC region completely. 



D. Effect of Vector Meson Coupling 

Now we ask the question whether there exists some 
mechanism that favors the diquark BEC. As we have 
seen, the diquark BEC happens in the beginning part 
of color superconductivity where the broken chiral sym- 
metry starts to restore and the effective quark mass is 
still large enough. As proposed in (65j , the quark inter- 
action in vector meson channel may be a candidate to 
slow down the chiral symmetry restoration and enhance 
the BEC formation. Now we include a new interaction 
term to the NJL lagrangian, 



£ v — G v 



(61) 



which corresponds to vector mesons. At finite density, a 
new condensate p v = 2Gv('^/ , 7o' ! / , ) which is proportional 
to the baryon number density should be considered. This 
condensate induces a new energy term — p v /4G v in the 
thermodynamic potential and it enters the gap equations 
by replacing /is/3 by p_b/3 — p v . The gap equation for 
p v at zero temperature reads 



d 3 p 

(2tt) 3 



5 

EZ 



e -e: 



(62) 



In Fig[l3] we calculate the effective quark mass and 
the diquark condensate at the coupling ratio r\ = 1 and 
for three values of the vector coupling G v . For G v = 
the diquark BEC region is almost not visible, but for a 
reasonably large vector coupling such as G v /G s = 0.3 
or 0.5, there appears a large diquark BEC region. It is 
clear that the vector meson channel really slows down the 
chiral symmetry restoration and favors the formation of 
diquark BEC. There may exist other mechanisms that 
favor BEC, such as the axial anomaly in QCD 66]. 
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FIG. 12: The diquark condensate A, effective quark mass m 
and the chemical potential as functions of baryon chemical 
potential for several values of the coupling ratio 77 at zero 
temperature. The color neutrality is considered. 



E. Chromomagnetic Instability 

For realistic QCD matter in nature such as in com- 
pact stars, beta equilibrium and charge neutrality should 
be considered. These effects induce a chemical potential 
mismatch 5/i between u and d quarks. Assuming the 
ground state to be a uniform phase, the mismatch ef- 
fect will lead to an interesting gapless color supercon- 
ductivity phases[67], [68j]. However, it was found that 
the gapless phase in the weak couplingregion suffers 
the so-called chromomagnetic instability |69|: The Meiss- 
ner masses squared of some gluons are negative. While 
the chromomagnetic instability indicates that the ground 
state of the superfluid may favor to be in some non- 
uniform phase such as LOFF phase [7(1 l7ll . the in- 
stability may be cured in some region. Based on a two 
species model 73] , it is found that the magnetic instabil- 
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FIG. 13: The diquark condensate A, effective quark mass m 
and the chemical potential /ijv as functions of baryon chemical 
potential for several values of the vector couplings G v and at 
a fixed diquark coupling ratio r\ — 1 at zero temperature. 



ity should be fully cured in the BEC region. However, 
the 4-7th gluons' instability is not yet examined. Here 
we try to complete this work. 

The analytical expressions of the Meissner masses 
squared TO4 for the 4- 7th gluons and m| for the 8th gluon 
are listed in Appendix [Bl We show the Meissner masses 
squared in FigO]|] in both BCS and BEC states. In the 
BCS case with m <C /z r , our result is consistent with 
the analytical expression 69] . When we approach to the 
BEC region, both the 4- 7th and 8th gluons' instabilities 
are partially cured. In the BEC region with m > fj, r , they 
are fully cured. Our investigation here is consistent with 
the study based on a non-relativistic model[74j]. Since 
the chromomagnetic instability is fully cured, the ground 
state in the BEC region should be a uniform phase. We 
should address that such a phenomenon is quite similar 
to what happens in non-relativistic systems [78l l79l. |8C| . 




FIG. 14: The Meissner masses squared, scaled by m 2 g — 
4a s /i^/(37r) with a a being the QCD gauge coupling constant, 
for the 4- 7th gluons (upper panel) and 8th gluon (lower panel) 
as functions of A/8[i. 



IV. PION BEC-BCS CROSSOVER AT FINITE 
ISOSPIN DENSITY 

Another BEC-BCS crossover in dense QCD happens at 
finite isospin density where pions become condensed 
At small isospin density, the QCD matter is a BEC of 
charged pions, but at ultra high isospin density, decon- 
finement happens and the matter turns to be a BCS su- 
perfluid with quark-antiquark cooper pairing^. There- 
fore, there should be a BEC to BCS crossover when the 
isospin chemical potential increases. From the similarity 
between pions in three color QCD and diquarks in two 
color QCD, the pion superfluidity discussed in this sec- 
tion is quite similar to the diquark superfluid in Section 
HIl We start from the two flavor Nambu-Jona-Lasinio 
model with only scalar meson channel, 

C = - m ) iJ; + G s [(^) 2 + (^75tV-) 2 ] • (63) 

The key quantity describing the system is the partition 
function 



where the isospin chemical potential /!/ corresponds to 
the third component / 3 of the isospin charge. Using the 
Hubbard-Stratonovich transformation we introduce the 
auxiliary meson fields a and 7T, and the partition function 
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can be written as 

Z = J [di>] [diP] [da] [dir]ef° dT 1 d3x£rff (65) 
with the effective Lagrangian 



4G., 



(66) 



where the kernel K, is defined as 

A*/ 

^ - "to -r — 

Integrating out the quark degrees of freedom, we obtain 



K[a, 7r] = ij^dn - m + — 70^3 - (<r + 1751" ■ tt)- (67) 
the quark degrees of 
Z = j[do\[dK)e- s ^ a '' (OS) 
effective ac 

S eff [<r,7r]= /"dr / rf 3 x " ^ " - Tr In K[a, tt], (69) 



with the meson effective action 



A. Pion Fluctuation at T > T c 

As we expected, when the isospin chemical potential 
fii becomes larger than the pion mass m ff in the vacuum, 
the charged pions will condense at low temperature. In 
this subsection we focus on the region above the criti- 
cal temperature T c where the pion condensate vanishes. 
After the field shift for cr, the effective action at zeroth 
order in meson fields gives the mean field thermodynamic 
potential Q, 



n = —s {0) = ( - m ~ TO °) 2 



1 



lndetS -1 , (70) 



/3V cS 4G S ' (3V 

where S is the quark propagator at mean field level 



The quadratic term of the effective action reads 



',< - m + y 70^3 



(71) 







dT 



d 3 x 



4G S 



+ -Tr{ l SS[a,7r] l SS[a, 7 r]}. 



(72) 



where S is defined as E[<7, tt] = a + ij 5 r ■ ir. Going to 
momentum space, the effective action can be evaluated 
as 



3(2) 
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2G 
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2Gl 
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2G~ S 



n ff (fc) 
n.o(fc) 



a{-k)cx{k) 
7r (-fc)7r (fc) 

7T + ( — fc)7T_(fe) 

7r_(-fc)7r + (fc), (73) 



where 7r± = (7Ti ± i-K 2 )/V^ are the positively and nega- 
tively charged pion fields, and the polarization functions 
for the mesons can be evaluated as 

n ff (fc) = N c a 1 (k;iJL I ), n T0 (fc) = A c n 2 (fc ;Ai/ ), 
n T+ (fc) = jv c n 3 (fc ;w ), n T _(fc) = iv c n 4 (fc ;M/ ) (74) 

where iV c is the color degree of freedom. 

Without loss of generality, we consider the case m > 
where the 7r + mesons become condensed at low temper- 
ature. The transition temperature T c is determined by 
the well-known Thouless criterion 



1-2G S IL K (fc = 0,k = 0) 







T=T C 



(75) 



together with the gap equation for the effective quark 
mass in derived from the first order derivative of the ther- 
modynamic potential with respect to m, 



m — mo 
8A C G, 



ni 



d 3 P l-f(E+)-f(E-] 



(2tt) 3 



(76) 



T=T C 



with Ep = E p ± fii/2. The critical temperature as a 
function of isospin chemical potential is exactly the same 
as in Fig[T] if we replace the baryon chemical potential 
l^B by the isospin chemical potential /i/. The diquark 
condensed phase starts at \xi = m ? , and the critical tem- 
perature can be well described by[60j 



T r =Tn 



Mi 



(77) 



where To is again the temperature of chiral symmetry 
restoration at /i/ = 0. 

The energy dispersions of the mesons are defined by 
the poles of their propagators, 



l-2GJT(fc = w(k),k) = 0, 



cr,7r ,7r + ,7r_. (78) 



Since the sigma and neutral pion do not carry isospin 
charge, they do not obtain an isospin chemical poten- 
tial and their masses are defined as m a — w CT (0) and 
m^a = U! V0 (Qi). However, for masses of charged pions, 
we should subtract the corresponding isospin chemical 
potential from the dispersions, namely we take m w+ = 
k-V + (0) + A*/ and m 7r _ — u lx _ (0) — /ij. For [i\ > 0, at the 
transition temperature T c , the 7r+ mass is equal to the 
isospin chemical potential, 



m 7r+ (T c ) = Hj. 



(79) 



The numerical results for meson masses and dissociation 
temperatures are the same as in Figs|3] and [3] if we take 
the correspondence d <-> 7r + , d <-> 7r_ , 7r <-> ttq and a a 
between the two cases. 

Similarly, we can investigate the spectral function 
p(u, k) for 7T + , 



p(w,k) = -2ImD R (a;,k), 



(80) 
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where Db(w, k) = T>(uj + irj, k) is the analytical continu- 
ation of the pion Green function T>(iLu ni k) defined as 



T>(iu) n ,k) 



2G S 



i-2G s n 7r+ (^„,k)- 



The spectral function is also similar to that of diquarks 
in Section |TJ For /i/ < 2m, ir + can decay into a quark- 
antiquark pair only when its energy satisfies lo > 2m— /i/, 
otherwise it remains a stable bound state. However, for 
fii > 2m, it can decay into a quark- antiquark pair at 
any energy and hence becomes an unstable resonance. 
The numerical result for the spectral function is similar 
to that shown in FigJU In conclusion, with increasing 
isospin chemical potential, the pions above the critical 
temperature become more and more unstable and finally 
disappear, which indicates a BEC-BCS crossover. 

We can also define an effective non-relativistic chemical 
potential fiN = Hi/2 — m - In the deep BCS limit with 
fij — ► oo, we have m — ► and /zjv — * In the deep 

BEC limit with [n — > m^, m tends to be its vacuum 
value, and [i^ becomes negative and its absolute value 
approaches to half of the pion binding energy 2m — m w . 
The effective quark mass m and the chemical potential 
/xjv at the critical temperature are the same as in Figj5l 
and the isospin chemical potential /ij for the crossover is 
still about 240 MeV. 



B. 



Pion Condensation at T < T c 



At low enough temperature, the pions become con- 
densed, we should introduce not only the chiral conden- 
sate (a) but also the pion condensates, 



A. 

7T 



V2 



A nonzero pion condensate A T ^ means spontaneous 
breaking of Us(l) symmetry corresponding to the isospin 
charge generator I3. The phase factor 9 indicates the 
direction of the U3(l) symmetry breaking. For a homo- 
geneous superfluid, we can choose 9 — without lose of 
generality. A gapless Goldstone boson will appear, which 
can be identified as the quantum fluctuation in the phase 
direction. 

After a field shift a — > (a) + a and 7Ti — > A^ + tti , the 
effective action at zeroth order in meson fields gives the 
mean field thermodynamic potential, 



n = 



1 



5 off 



(m - m ) 2 + A 2 



1 

w 



lndetS^ 1 , (83) 



where is the quark propagator in the pion superfluid 
at mean field level and can be evaluated as a matrix in 
flavor space, 



Suuip) S ud {p) 
S du (p) Sdd(p) 



(84) 



with the elements 

(iv n + E~ ) A+7o (iv n - E+) A_ 7o 



Suu 
Sdd ■ 
S u d 
Sdu 



[iv n - E~) A_7o {iv n + E+) A + 7o 



{i Vn y-{E^y 

-iA w A+75 



-iA 7r A_7 5 



-zA 7r A_7 5 -zA^A+75 
{iv n f-(E-y (iv n ) 2 - {E+) 2 ' 



(85) 



where E^ = J (Ep) 2 + A 2 are quark energies. 

The momentum distributions of quarks and anti- 
quarks can be calculated from the positive and negative 



energy components of the diagonal propagators S u 

Sdd, 



and 



n u {p) 
n d (p) 



n d(p) = \ 



n u{p) = 2 



E n 
El 



!--T 

Ett 



• (86) 



The gap equations to determine the effective quark 
mass m and pion condensate A^ can be obtained by min- 
imizing the thermodynamic potential, 



on 

dm 



0. 



on 



0. 



(87) 



With the explicit form of the thermodynamic potential 



n = 



(m~m ) 2 + Al 
4G S 



-4A^ r 



d 3 p 

(2Vp 



[aE+) + aE-)] 



(82) we obtain the explicit gap equations at zero temperature, 



m — mo 

A-n- 



AN c G s m 



4N c G s Av 



d-'> 



P 1 



-1^ 

(2^) 3 £p \EZ 

d 3 p 



4 

E-n 



(27T) 



11, 

W + F+ ] - (89) 



The numerical results for m and A T are the same as in 
FigEJ if we replace ug by \xi. The results agree well with 
lattice dataflfl EH G2] at least at small isospin chemical 
potential. For /ij < m v , the ground state is the same 
as the vacuum state and the isospin density keeps zero, 
while for fij > m v , the pion condensate and isospin den- 
sity become nonzero. In the isospin chemical potential 
region above but close to the critical value fii = m w , 
the effective quark mass and pion condensate can be well 
described by 



m(0) 

A(/zi) 
M(0) 



W(0) 



= \ 1- 



Mi 



(90) 
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The explanation of the BEC-BCS crossover at zero 
temperature is similar to the diquark case in Section HT1 
For < fJ-i/2 < m(fj,i), the minimum of the dispersion 
E~ is at |pj = where the energy gap is \/ [i z N + Aj. 
For large enough fij, the minimum of the dispersion is 
shifted to |p| ~ M//2 where the energy gap is A. The 
momentum distribution for u and anti-d quarks is the 
same as in FigJTl which shows a BEC-BCS crossover at 
finite isospin chemical potential. Similar to the diquark 
case, the chemical potential /Jj for the crossover is 



[2m(0)m»] 



(91) 



which is about 230 — 270 MeV when the parameter values 
change reasonably. 

The second order effective action which is quadratic in 
the meson fields and controls the meson behavior can be 
evaluated as 



dr / fx 

o J 4G S 

^Tr{S*E[*,ir]S*H[tT,ir]}. (92) 



In momentum space, it is related to the meson polariza- 
tions, 



2G. S 



Ikj(k) 



M-QMQ (93) 



with i,j = a, 7r+ , 7T_ , 7To , where the polarization functions 
d 4 p 



liij are defined as 



Hij(k) = iN c 



(20' 



-Tr[T a S„{p + k)T b S„{p)]. (94) 



The masses of those new eigen modes in the superfluid 
phase are determined by the pole equation 



det 



^--ny(fc = M,k=0) 



= 0. 



(95) 



It can be shown that the neutral pion does not mix with 
the other mesons and is still an eigen mode of the system, 
but sigma and charged pions are strongly mixed, and it 
is the mixing that results in a gapless Goldstone boson. 
From the proportional relations H a -K + > IFj.^ ~ mA n and 
IFr + u-_ ~ A^[30j, the mixing between sigma and charged 
pions is very strong in the BEC region where m and A v 
are both large and coexist, but can be neglected at large 
isospin chemical potential. 

In summary, the phase diagram of QCD at finite 
isospin density is shown in Fig |15l At low temperature 
and low isospin chemical potential, the matter is in nor- 
mal hadron gas with chiral symmetry breaking. With 
increasing temperature, there should be a phase tran- 
sition from hadron gas to quark gas, indicated by the 
thick dashed line. When the isospin chemical potential 
becomes larger than the pion mass in the vacuum, the 
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FIG. 15: The proposed phase diagram of pion condensation 
in the T — Hi plane. 



pion BEC appears and keeps until another critical isospin 
chemical potential which is indicated by the verti- 
cal dashed line. At high enough isospin chemical poten- 
tial, the matter will enter the BCS superfluid where the 
quark-antiquark Cooper pairs are condensed. Between 
the BEC and BCS, there should exist a large crossover 
region. Since QCD at finite isospin chemical potential 
can be successfully simulated on lattice, such a BEC- 
BCS crossover can be confirm by measuring the quark en- 
ergy gap and comparing it with the pion condensate, and 
the pseudogap phase can be confirmed by investigating 
the quark spectral function. When finite baryon density, 
charge neutrality and weak equilibrium are taken into ac- 
count, they may have significant effect on the BEC-BCS 
crossover of pion condensation^. l76l |77| . 



V. CONCLUSIONS 

We have studied the BEC-BCS crossover in QCD at fi- 
nite baryon or isospin density in the NJL model at quark 
level. We investigated the BEC-BCS crossover in two 
aspects: (1) Above the critical temperature of the su- 
perfluid, diquarks or mesons are stable bound states at 
low chemical potential but become unstable resonances 
at high chemical potential; (2) At zero temperature, the 
effective non-relativistic chemical potential, dispersions 
of fermion excitations and the fermion momentum dis- 
tribution behavior significantly differently in the low and 
high chemical potential regions. The diquark BEC-BCS 
crossover in two color QCD at finite baryon density and 
the pion BEC-BCS crossover in real QCD at finite isospin 
density can be identified, since the quark confinement is 
not important in two color case. We expect that such 
a crossover can be confirmed in the lattice simulations, 
since the results from the NJL model agree quite well 
with the lattice data obtained so far. However, for real 
QCD at finite baryon density, whether there exists a di- 
quark BEC-BCS crossover is still an open question, since 
the confinement in this case is quite important. 

An interesting and important phenomenon we found in 
this paper is that the BEC-BCS crossover we discussed 
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is not induced by simply increasing the coupling con- 
stant of the attractive interaction but by changing the 
corresponding charge number. During the change of the 
number density, the chiral symmetry restoration plays an 
important role in the study of BEC-BCS crossover. 

When an chemical potential or fermion density mis- 
match between the two pairing species is turned on, the 
BEC-BCS crossover will be dramatically changed. The 



system may go from some non-uniform phases such as 
LOFF in the BCS reg ion to some unifor m g apless phase 
in the BEC region^ O US EE E2, ES E3|- 
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APPENDIX A: POLARIZATION FUNCTIONS ni,n 2 ,n 3 ,n 4 



The meson and diquark polarization functions ETi, Ila, II3 and n4 above the critical temperature for diquark or pion 
condensation can be evaluated as explicit functions of temperature T and corresponding chemical potential fx, 



n x (fc; M ) = - 



n 2 (fc ;M ) = - 



d 3 p 

(2tt) 3 



d 3 p 

(2tt) 3 



TZ 



fp fp /q /q I p ~t~ fp fq /q 

fc - £ q + E p k + Eq - E p 

o _ f- — f+ — f- — f+ 9 — f~ — f+ — f~ — f+ 

Jp Jp Jq Jq Jp Jp Jq Jq 



ko — Eq — Ep ko + Eq + E p 

3 p / p /q f C{ f p f p I q ^q \ rj — h 

ko-E cl + Ep fc + E^ - E p ) 

_ /- _ f+ _ f- _ f+ 2 — f~ — f+ — f ~ — f + 

Jp Jp Jq Jq Jp Jp jq Jq 



TZ 



II 4 (fc;/i) = -2 



(2tt) 3 



d 3 p 



ko — Eq — Ep 

f+ _ f- 

J p ./ q 



k + Eq + Ep 

f- _ f+ 

J p J q 



fc + (i - Eq + E p ko + (i + Eq - E p 

1 - fp - /q 1 - /p + - /q + 



k + (X - Eq - Ep k + fi + Eq + Ep 



fp i q 



f+ _ f- 

J p J q 



ko — f-i — Eq + Ep ko — fi + Eq — Ep 



+ 



^q -r ^p 

+ - f+ 



1 - f; - / q 



1 - fp - /q 



k - fi - Eq - E p k - [i + Eq + Ep 



TZZ 



T 



TZ 



(Al) 



where we have defined q = p+k, E p = \/p 2 + m 2 , fp = f(E p ±fi/2) with f{x) = (l + e x ^ T ^j 1 being the Fermi-Dirac 
distribution function and T^ = 1 ± (p • q =p m 2 ) / (E p Eq). 



APPENDIX B: MEISSNER MASSES SQUARED OF GLUONS 



In this Appendix we list the explicit form of the Meissner masses squared to 2 for the 4- 7th gluons and m\ for the 
8th gluon in two- flavor color superconductivity phase which are quoted from (85l . Ell > 



2g 2 r d 3 p 
~J (2^)3 
2.9 2 f d 3 p 



3-^)A(p) + ^B(p) 



E l 



C{P) 



E l 



D{P) 



(Bl) 
(B2) 
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with the QCD gauge coupling constant g and the functions 

= (E A ) 2 -E p E+ + A 2 e(- El ) + e(-E 2 )-l (E+) 2 -E p E+ + A 2 ®(-E 3 ) + e(-E i )-l 1 
[P> {E- A Y-{ElY E A {E- A Y~{ElY El E p 

B(p) = -8(Eh.) - S(E 2 ) - S(E 3 ) - 6(Ei), 

_ 2 f e(-E 5 )-e(E 2 ) e(-E 7 ) - ejE^ 2 { q(-e 5 ) - e^) , g(-e?) - q{-e 2 ) 

<^(P) - u - \ ^ ; ^ 1 B i S ) +w 

D(p) = u 2 _ 



E5 + E 2 E-j + Ex J \ E5 — Ei Ej — Ei 

e(-E 6 )-Q(E 3 ) , Q(-E 8 )-e(E 4 )\ , a /6(-£? 6 ) - 6(-E 4 ) , 9(-£ 8 ) - Q(-E 3 ) \ 2 



+ E 3 Eg + Ei J \ E 6 — Ei E s — E 3 J E p 

&{-e 6 ) - e(-E 2 ) , &(-e 8 ) - &(-Ei)\ , 2 /ec-^-e^) , 9(-£ 8 ) - 9(£ 2 y 



i?6 — -E2 £>8 — Ei J V £>6 + -Ei -£"8 + -E; 



^2 

2 / 9(-g 5 )-e(-g 3 ) 9(-£ 7 )-9(-£ 4 ) \ 2 / 9(-£ 5 )-9(e 4 ) e(-£ 7 )-e(£ 3 ) \ , 
+ V s 5 -b 3 £7 -Ei ) +\ e 5 + e 4 e- + e 3 )' [ ' 

where the quark energies are defined as 

Ei = E A - Sfi, E 2 = E A + <5/i, E 3 = E+- Sfi, E 4 =E++ Sfi, 

E 5 = E+-Sfx, E 6 = E^+8fJi, E 7 = E+ + Sfx, Eg = E^ - Sfi (B4) 

with the chemical potential difference S/i between u and d quarks and E A and E^ being listed in Section IIII1 the 
coherent coefficients u± and v\ are defined as u± = (l + E±/E%) /2 and v 2 ± = (l - E^/E^) /2. Note that we have 
added the terms 1/E P to A and 2/E p to C to cancel the vacuum contribution |73| . In this way the Meissner masses 
squared are guaranteed to be zero in the normal phase with A = 0. 
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